Effect of orbital-overlap dependence in density functionals 
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The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation 
functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies 
simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy 
density as an input. We study for the first time the effect and importance of the dependence of 
MGGA on the kinetic energy density through the dimensionless inhomogeneity parameter, a, that 
characterizes the extent of orbital overlap. This leads to a simple and wholly new MGGA exchange 
functional, which interpolates between the single-orbital regime, where a = 0, and the slowly varying 
density regime, where a ~ 1, and then extrapolates to a -^ oo. When combined with a variant of 
the Perdew-Burke-Erzerhof (PBE) GGA correlation, the resulting MGGA performs equally well for 
atoms, molecules, surfaces, and solids. 
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Kohn-Sham (KS) density functional theory (DFT) 
[l|,[3| is one of the most widely used methods in condensed 
matter physics and quantum chemistry. In this theory, 
the exchange-correlation energy E^c as a functional of the 
electron spin densities Ti-|-(r) and ?T-4,(r) must be approx- 
imated. Semilocal approximations (e.g., Refs. |3l-[l0|) of 
the form 

^xc[nt,n-i]= d^rne^cin^,ni,\7n^,\/ni,T^,Ti) (1) 

require only a single integral over real space and so 
are practical even for large molecules or unit cells. In 
Eq. ([1]), Vn-|-,4, are the local gradients of the spin densi- 
ties, T^^i = ^^ |VV'/ct,4.| /2 the kinetic energy densities 
of the occupied KS orbitals il^ka- of spin a, and e^c the ap- 
proximate exchange-correlation energy per electron. All 
equations are in atomic units. Semilocal approximations 
can be reasonably accurate for the near-equilibrium and 
compressed ground-state properties of ordinary matter, 
where neither strong correlation nor long-range van der 
Waals interaction is important. They can also serve as a 
base for the computationally more-expensive fully nonlo- 
cal approximations needed to describe strongly correlated 
systems [11| and soft matter [l2| . 

The meta generalized gradient approximation 
(MGGA) is the highest semilocal rung of the so-called 
Jacob's ladder in DFT [Ij]. In addition to the spin 
densities n-|- and n^ that are used in local spin density 
approximation (LSDA) [l|, 0, |j] and their gradients 
that are further included in generalized gradient ap- 
proximation (GGA) 0-01, MGGA includes the kinetic 
energy density that can be used, as in the revised 
Tao-Perdew-Staroverov-Scuseria (revTPSS) MGGA fl, 
to distinguish the single-orbital regions from the orbital- 
overlap regions. However, the dependence of MGGAs 
on the kinetic energy density is understood much less 
than that on the density gradient, which MGGA inherits 
from GGA. Such understanding is highly demanded by 
the construction of not only MGGAs theirselves but 
also fully nonlocal approximations that are based on 
MGGAs. Therefore, this could be largely beneficial to 



expediting the shift in DFT from the dominant GGAs 
to the generally more accurate and computationally 
comparable MGGAs @, [iJI- In this article, we will show 
the importance and effect of the r-dependence, leading 
to a new MGGA that respects a tight Lieb-Oxford 
bound [l^, [la| and performs equally well for atoms, 
molecules, surfaces, and solids. 

The semilocal exchange energy of a spin-unpolarized 
density can be written as: 



El 



= d 



\n)F^{p,a). 



(2) 



Here, e™'^(ri) = —-^{^Y^^/fs is the exchange en- 
ergy per particle of a uniform electron gas with r^ = 
(47rn/3)-i/3, p = |V7i|V[4(37r2)2/3„8/3] ^ ^2^ ,^^^ 

a = {t - t^)/t'''^^^ ^ 0. In the latter, r = Y.a'^'^^ 
T^ — i I Vnp/n is the von Weizsacker kinetic energy den- 
sity, and r""'^ — r^(37r^)^/'^n^/'^ is the orbital kinetic en- 
ergy density of the uniform electron gas. The expression 
for the s pin- polarized case follows from the spin-scaling 
relation [l3] . The enhancement factors F^ (which equals 
1 in LSDA) of GGAs, independent of a, are often made 
monotonically increasing with s as in the standard PBE 
GGA |5| , or for a large range of s [l8l42Cl| , and therefore 
favor less compact systems than LSDA does (e.g., lower- 
ing the energy of a collection of dissociated atoms relative 
to that of the molecules, or lowering a surface energy, 
or enlarging lattice constants of solids). The revTPSS 
MGGA includes the extra ingredient, a, and recovers 
the exact exchange energy of the ground state density 
of the hydrogen atom and the finiteness of the exchange 
potential at nuclei, where a = (single-orbital regime). 
Then, at a w 1 (slowly- varying density regime), it re- 
stores the second order gradient expansion for a wide 
range of density and further recovers the fourth order 
gradient expansion of a slowly varying density. There- 
fore, we believe that the revTPSS enhancement factor 
is accurate for small s around a k 1. However, in the 
construction of revTPSS, there is no other constraint to 
guide the functional approaching from a = to a w 1. 
revTPSS also has an order-of-limits problem [2l|, [22] — 



the enhancement factor has different values when differ- 
ent sequences of the Hmits s — > and a — >■ are taken, 
as shown in Fig. 



1(a) 



Here, we propose a simple exchange enhancement fac- 
tor that disentangles a and p by a means of separability 
assumption, 

F^ (P, «) = Fl{v) + /(«)[^°(P) - Fl{v)\ (3) 

where Fl{v) = Ff\'p, a = 1) = 1 + k - k/(1 + fi'^^p/n) 
andi^O(p) = F^"'(p,a = 0) = l-hK-Av/(l-F(At°'^p-hc)/K). 
F"'' {p, a) interpolates between F° (p) and F^ (p) through 
/(a) = (1 — a^)^/(l + 0"^ +ck^), which is chosen to guar- 
antee for the functional the second order gradient ex- 
pansion, good exchange jellium surface energies, and the 
Hartree-Fock exchange energy of the 12-electron hydro- 
genic density (ls^2s^2p^3s^) with the nuclear charge Z = 
1 [23] {Z — Im the following discussions if not mentioned 
otherwise) . See Ref. [2J for the derivatives of F™^ (p, a) 
with respect to p and a which are needed for the self- 
consistent implementation. For a slowly varying density 
where a « 1, the second term of the left hand side of Eq.[2] 
is negligible and of order V^ . F^"* {p, a) then reduces to 
F,}{p) and recovers the second order gradient expansion 
with the first principle coefficient /i*^^ = 10/81 f25'| as 
in PBEsol ,6j. c = 0.28771 and n = 0.29 are two pa- 
rameters fixed by the exchange energies of the hydrogen 
atom where a = and the 12-electron hydrogenic density 
which is used to guide the functional from a — Oto a ^ 1. 
They also deliver excellent exchange energies for other 
hydrogenic densities (see Ref. |24I ). No point in the (c, 
k) parameter space could be found without violating the 
loose Lieb-Oxford bound [H (ie., k sC k^^ = 0.804 0) 
if we use a spin-unpolarized hydrogenic density with elec- 
tron number less than 12. The obtained k — 0.29 sug- 
gests a tight Lieb-Oxford bound and leads to a very 
fiat F^{p) as shown in Fig. 1(b) The resulting small 

derivative at nuclei, — ^ d^'"~ ls=o.376 — 0.22, comes 
close to satisfying the finiteness constraint on the ex- 
change potential at nuclei. Compared to revTPSS, this 
much simpler form doesn't have the order-of-limits prob- 
lem while recovering regions of small s around a ~ I oi 
revTPSS as shown in Fig. 1(a) In the following discus- 



sion, we combine this exchange functional with the vari- 
ant of the PBE correlation (denoted as vPBEc), where 
/3(rs) = 0.066725(1 + 0.1r^)/(l 4- 0.1778rs) as used in 
revTPSS Q. Although the vPBEc is not one-electron 
self-correlation free, its error is small (about 0.006 Ha for 
the H atom), which is usually largely cancelled out of 
atomization energies involving H. It has an accurate cor- 
relation energy for a two-electron ion of nuclear charge 
Z ^ oo, which is -0.0479 Q, better than -0.0510 from 
TPSS and -0.0527 from revTPSS in comparison with the 
exact value -0.0467 [2J|. And it also has accurate corre- 
lation jellium surface energies Q. The resulting MGGA 
respects a tight Lieb-Oxford bound with Fxc ^ 2.137, 
while the loose Lieb-Oxford bound is 2.273 [il,[ii|. 

From the experience of GGAs [a,[a], we know the faster 
an enhancement factor grows with s, the more preference 





(b) 



FIG. 1: Exchange enhancement factors vs. a (a), and s (b). 



of the functional towards less compact systems. Since 
the enhancement factor of the present exchange func- 
tional as shown in Fig. 1(b) is largely depressed for a 



large range of s compared to that of revTPSS, it then 
seems to be a reasonable guess from the experience that 
the present exchange functional with the vPBEc would 
give too small lattice constants for solids and too high 
atomization energies for molecules. However, our results 
(given later) show a different scenario that our MGGA 
performs equally well for atoms, molecules, surfaces, and 
solids. The seeming contradiction between the prefor- 
mance and Fig. 1(b) is well resolved by resorting to the 
a-dependence shown in Fig. l(a)[ which is much less un- 
derstood. In previous constructions and analysis of MG- 
GAs [i, H, [Tj], it is emphasized that a enables MGGAs 
to distinguish the single-orbital regions from the orbital- 
overlap regions. Here, we further stress that the mono- 
tonically decreasing dependence of an enhancement fac- 
tor on a is qualitatively equivalent to the monotonically 
increasing s-dependence, which explains the seeming con- 
tradiction and will be rationalized by the following two 
observations. 

The first is on the changes in s and a distributions 
from the 10-electron hydrogenic density {ls^2s^2p^) to 
the 12-electron one {ls^2s^2p^3s^) ^, and on the cor- 
relation between the changes. Note these two densities 
are spherically symmetric. Fig. [5] shows the s and a 
distributions of these two densities. The shell structure 
can be easily recognized and roughly identified by a with 
a < 1 for shell regions and a > 1 for intershell regions. 
We can tolerate the confusion caused by this definition 
for the tail regions, where a could be or ^ 1 as shown 
in Fig. [21 since the tail regions are energetically less im- 
portant. When an electron is in the intershell region, and 




FIG. 2: The s and a distributions for the 10- and 12-electron 
hydrogenic densities. The line y=l is plotted to help roughly 
identify the shell (a < 1) and intershell (q > 1) regions. 



a is large even where s is sniall, the electron's exchange 
hole is probably not centered close to the electron, but 
is spread out over the smaller inner shell and the larger 
outer shell. This spreading of the hole will make the 
exchange energy density in the intershell region less neg- 
ative than it would be for a slowly-varying or uniform 
density. One can imagine i^^; < 1, as can happen in the 
present exchange functional for small s and large a that 



is shown in Fig. 1(a) 



When two 3s hydrogenic electrons are added to the 10- 
electron hydrogenic density, part of the outermost shell 
region (a < 1) changes into the intershell region (a > 1), 
which is associated with a decrease of s. Although it's 
not certain that an increase of a is always associated 
with a decrease of s during the formation of the inter- 
shell of an atom, it's very likely for the intershell region 
between the outermost core and the valence of an atom 
within a solid, an important region for determining the 
lattice constants of solids [2y, [27| (as will be discussed in 
the second observation). Then, the presumable correla- 
tion between s and a in the intershell regions of a solid 
suggests that monotonically decreasing a-dependence of 
an enhancement factor has the qualitatively same effect 
as monotonically increasing s-dependence does for these 
regions. 



TABLE L Error statistics of lattice constants (A) of the SL20 
solids J28| ] and atomization energies (kcal/mol) of the AE6 set 



See the text for the definitions of F^ 



F* 



and F" 





LSDA PBEsol 


^Omt 


^*mt 


F^' 


ME 
SL20 MAE 


-0.081 -0.012 
0.081 0.036 


-0.079 
0.079 


-0.126 
0.126 


0.016 
0.023 


ME 
AE6 MAE 


77.4 35.9 
77.4 35.9 


55.1 
55.1 


22.0 
22.6 


0.6 
5.5 



The second observation concerns the variations of the 
lattice constants of the set of 20 solids (SL20) [2^ and 
the atomization energies of the AE6 molecule set [2^ 
in response to changes of the a dependence of the en- 
hancement factor. In order to show the dramatic effect of 



the a dependence on the structural and thermochemichal 
properties and to deduce its origin, f{a) is compared 
to two variants. The first variant is /°(a) — 0, which 
is independent of a and actually results in the PBEsol 
GGA but with different k and (3{rs)- The second one is 
f*{a) = /(a) * hia) with h{a) = 2/[e("-i)/'' + 1] - 1. 
h(a) is equal to 1 for a < 1 and to -1 for a > 1, re- 
spectively, if 7 ^- 0. Here, we choose 7 = 0.1 for nu- 
merical reasons. Compared to /(a), /*(«) flips the a 
dependence for a > 1 from monotonically decreasing to 
monotonically increasing, and therefore favors more the 
regions with a > 1. The choice of the demarcation point 
at a = 1 is natural in view of the first observation and 
because it helps to satisfy the second order gradient ex- 
pansion, /(a) and the two variants — whose curves as 
functions of a are given in Ref. |2J — result in three dif- 
ferent enhancement factors, F™*, ^i?'"', and f^'"*. 



TABLE II: Error statistics for various density functionals. 
The LSDA and PBE enthalpies of formation are from Ref. |y. 
The lattice constants for M06L are from Ref. 30i, where 18 
solids, mostly overlapping with the SL20 set, were tested, a: 
Ref. [33 



LSDA PBE 



M06L 



revTPSS present 



Exchange energies (Ha) of rare gas atoms 
ME 2.274 0.219 0.204 0.291 0.068 

MAE 2.274 0.219 0.210 0.293 0.111 

Atomization energies (kcal/mol) of the AE6 molecules 
ME 77.4 12.4 3.2 3.3 0.6 

MAE 77.4 15.5 4.2 5.9 5.5 

Dissociation energies (kcal/mol) of the W6 water clusters 
ME 5.2 0.0 -0.2 -1.0 0.0 

MAE 5.2 0.3 0.4 1.0 0.1 

Enthalpies of formation (kcal/mol) of the G3 molecules 
ME -121.9 -21.7 -1.6 -3.6 -1.6 

MAE 121.9 22.2 5.2 4.8 8.3 

Jellium surface exchange energies (%) 
MRE 45.8 -20.9 -75.9 -1.0 -7.3 

MARE 45.8 20.9 75.9 2.2 8.0 

Jellium surface exchange-correlation energies (%) 
MRE -0.4 -3.1 24.5 2.6 -0.3 

MARE 0.4 3.1 24.5 2.6 1.6 

Lattice constants (A) of the SL20 solids 
ME -0.081 0.051 0.015" 0.015 0.016 

MAE 0.081 0.059 0.071" 0.033 0.023 



TableUshows the mean error (ME) and the mean abso- 
lute error (MAE) of the lattice constants of the SL20 set 
[28| . and the atomization energies of the AE6 set [20| 
from LSDA, PBEsol, and variants of F™'. The alle- 
viation, from LSDA to F°'"^^ and then to F^'^', of the 
overestimation in the atomization energies and of the 
underestimation in the lattice constants, suggests that 
the built-in monotonically increasing s-dependence and 
monotonically decreasing a-dependence in the enhance- 
ment factors reduce the preference of LSDA towards com- 
pact systems. However, from F"* to F^'"', where only 



the monotonicity of the a-dependence in the range of [1, 
oo] is flipped, the sohds in the SL20 set are drastically 
shrunk to such a surprising degree that the lattice con- 
stants of i^^'"' are significantly smaller than even those 
of LSDA. Since the important region in terms of deter- 
mining the lattice constant of a solid for a functional has 
been identified [2a, [23] to be the intershell region of the 
constituent atoms between the outermost core and the 
valence regions, the drastic shrinkage from F^"' to F*'"* 
is a strong indication that the sheh and intersheU re- 
gions are associated with q < 1 and a > 1, respectively. 
Compressing a sohd turns part of the outermost core 
and the valence regions (a < 1) into intershell regions 
(a > 1) between them, which F*'"' favors more than 
F^"' does. The absence of the monotonicahy decreasing 
a-dependence in F^?'"*, leading to the shrinkage of the 
solids, could be compensated by enhancing the mono- 
tonically increasing dependence on s, as PBEsol does by 
using K^'^ . This implies a decrease of s during the for- 
mation of the intershell regions and thus corroborates the 
correlation between s and a during the formation of the 
intershell region observed in Fig. [51 Therefore, in F"*, 
both the monotonically increasing s-dependence and the 
monotonically decreasing a-dependence have the effect 
of penalizing the formation of the intershell regions and 
enlarging the lattice constants. Similar deterioration is 
also found for the atomization energies of the AEG set for 
F*''^* compared to F^"*. Unlike for solids, F*'"' still sig- 
nificantly improves the atomization energies of the AE6 
set over F"'"', and thus LSDA, suggesting that the a- 
dependence in the range of [0, 1] has stronger influence 
in atoms and molecules than in solids. Remarkably, the 
monotonically decreasing a-dependence used in F"* sig- 
nificantly improves the overestimated atomization ener- 
gies of the AEG set of F°'"' to an excellent accuracy level 
with MAE of 5.4 kcal/mol. This implies that monotoni- 
cally decreasing a-dependence is in general able to make 
a functional favor less compact systems, as does mono- 
tonically increasing s-dependence. 

Now, let's turn to the results for atoms, molecules, sur- 
faces, and solids, which are summarized briefly in Table [U 
in terms of the ME and MAE, or their relative analogs 
MRE and MARE. See Ref. M for full details of Table [fil 
Table [ni shows that the use of the exact exchange energy 
of the 12-electron hydrogenic density guarantees excellent 
exchange energies for atoms, resulting in good atomiza- 
tion energies for this simple functional. In all categories 
shown in Table [Til the present functional outperforms 
the standard PBE GGA (PBE is shghtly better than 
the present functional for the cohesive energies of the 
20 SL20 sohds as shown in Ref.'2i). Within the MGGA 



level. The heavily parameterized MOGL [T3| predicts ex- 
cellent atomization energies and dissociation energies for 
the W6 set, at which it aims during the construction. 
However, it is significantly wrong for the jellium sur- 
face energies contributed from the exchange and correla- 
tion terms, separately or together. The too-large MOGL 
jellium surface exchange-correlation energies imply that 
metal bulks are overstabilized, consistent with the too- 
small lattice constants of main group simple metals, e.g., 
Na and Al \30\. The present functional and revTPSS are 
more balanced for different categories and therefore more 
robust. Compared to revTPSS, the present functional is 
better for the SL20 solids, comparable for the jellium 
surface exchange-correlation energies but worse for the 
exchange part alone, and worse for the G3 molecules. 
However, the present functional predicts the most accu- 
rate dissociation energies of the W6 water clusters among 
the functionals, implying a good description for the hy- 
drogen bond. 

In summary, we have for the first time stud- 
ied the effect and importance of the dependence of 
computationally-efficient semilocal MGGAs on the ki- 
netic energy density through the dimensionless inhomo- 
geneity parameter a, and presented a new MGGA ex- 
change functional that disentangles a from the reduced 
density gradient s by the means of separability assump- 
tion. By varying the a— dependence in the exchange 
functional, we showed that the formation of the inter- 
shell region between the outermost core and the valence 
of an atom within a solid is associated with an increase 
of a and a decrease of s, suggesting that monotoni- 
cally decreasing a-dependence of an enhancement fac- 
tor is qualitatively equivalent to monotonically increas- 
ing s-dependence for these intershell regions. This has 
a significant impact on the construction of MGGAs and 
the MGGA-based nonlocal approximations, as exempli- 
fied by the present MGGA — which is overall comparable 
in performance, but quite different and much simpler in 
form, compared to the sophisticated revTPSS MGGA, 
and thus demonstrates the fiexibility and the rich struc- 
ture of MGGA brought by the extra ingredient of the 
kinetic energy density. 
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The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation 
functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies 
simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy 
density as an input. We study for the first time the effect and importance of the dependence of 
MGGA on the kinetic energy density through the dimensionless inhomogeneity parameter, a, that 
characterizes the extent of orbital overlap. This leads to a simple and wholly new MGGA exchange 
functional, which interpolates between the single-orbital regime, where a = 0, and the slowly varying 
density regime, where a ~ 1, and then extrapolates to a -^ oo. When combined with a variant of 
the Perdew-Burke-Erzerhof (PBE) GGA correlation, the resulting MGGA performs equally well for 
atoms, molecules, surfaces, and solids. 
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Kohn-Sham (KS) density functional theory (DFT) 
[l|,[3| is one of the most widely used methods in condensed 
matter physics and quantum chemistry. In this theory, 
the exchange-correlation energy E^c as a functional of the 
electron spin densities Ti-|-(r) and ?T-4,(r) must be approx- 
imated. Semilocal approximations (e.g., Refs. |3l-[l0|) of 
the form 

^xc[nt,n-i]= d^rne^cin^,ni,\7n^,\/ni,T^,Ti) (1) 

require only a single integral over real space and so 
are practical even for large molecules or unit cells. In 
Eq. ([1]), Vn-|-,4, are the local gradients of the spin densi- 
ties, T^^i = ^^ |VV'/ct,4.| /2 the kinetic energy densities 
of the occupied KS orbitals il^ka- of spin a, and e^c the ap- 
proximate exchange-correlation energy per electron. All 
equations are in atomic units. Semilocal approximations 
can be reasonably accurate for the near-equilibrium and 
compressed ground-state properties of ordinary matter, 
where neither strong correlation nor long-range van der 
Waals interaction is important. They can also serve as a 
base for the computationally more-expensive fully nonlo- 
cal approximations needed to describe strongly correlated 
systems [11| and soft matter [l2| . 

The meta generalized gradient approximation 
(MGGA) is the highest semilocal rung of the so-called 
Jacob's ladder in DFT [Ij]. In addition to the spin 
densities n-|- and n^ that are used in local spin density 
approximation (LSDA) [l|, 0, |j] and their gradients 
that are further included in generalized gradient ap- 
proximation (GGA) 0-01, MGGA includes the kinetic 
energy density that can be used, as in the revised 
Tao-Perdew-Staroverov-Scuseria (revTPSS) MGGA fl, 
to distinguish the single-orbital regions from the orbital- 
overlap regions. However, the dependence of MGGAs 
on the kinetic energy density is understood much less 
than that on the density gradient, which MGGA inherits 
from GGA. Such understanding is highly demanded by 
the construction of not only MGGAs theirselves but 
also fully nonlocal approximations that are based on 
MGGAs. Therefore, this could be largely beneficial to 



expediting the shift in DFT from the dominant GGAs 
to the generally more accurate and computationally 
comparable MGGAs @, [iJI- In this article, we will show 
the importance and effect of the r-dependence, leading 
to a new MGGA that respects a tight Lieb-Oxford 
bound [l^, [la| and performs equally well for atoms, 
molecules, surfaces, and solids. 

The semilocal exchange energy of a spin-unpolarized 
density can be written as: 



El 



= d 



\n)F^{p,a). 



(2) 



Here, e™'^(ri) = —-^{^Y^^/fs is the exchange en- 
ergy per particle of a uniform electron gas with r^ = 
(47rn/3)-i/3, p = |V7i|V[4(37r2)2/3„8/3] ^ ^2^ ,^^^ 

a = {t - t^)/t'''^^^ ^ 0. In the latter, r = Y.a'^'^^ 
T^ — i I Vnp/n is the von Weizsacker kinetic energy den- 
sity, and r""'^ — r^(37r^)^/'^n^/'^ is the orbital kinetic en- 
ergy density of the uniform electron gas. The expression 
for the s pin- polarized case follows from the spin-scaling 
relation [l3] . The enhancement factors F^ (which equals 
1 in LSDA) of GGAs, independent of a, are often made 
monotonically increasing with s as in the standard PBE 
GGA |5| , or for a large range of s [l8l42Cl| , and therefore 
favor less compact systems than LSDA does (e.g., lower- 
ing the energy of a collection of dissociated atoms relative 
to that of the molecules, or lowering a surface energy, 
or enlarging lattice constants of solids). The revTPSS 
MGGA includes the extra ingredient, a, and recovers 
the exact exchange energy of the ground state density 
of the hydrogen atom and the finiteness of the exchange 
potential at nuclei, where a = (single-orbital regime). 
Then, at a w 1 (slowly- varying density regime), it re- 
stores the second order gradient expansion for a wide 
range of density and further recovers the fourth order 
gradient expansion of a slowly varying density. There- 
fore, we believe that the revTPSS enhancement factor 
is accurate for small s around a k 1. However, in the 
construction of revTPSS, there is no other constraint to 
guide the functional approaching from a = to a w 1. 
revTPSS also has an order-of-limits problem [2l|, [22] — 



the enhancement factor has different values when differ- 
ent sequences of the Hmits s — > and a — >■ are taken, 
as shown in Fig. 



1(a) 



Here, we propose a simple exchange enhancement fac- 
tor that disentangles a and p by a means of separability 
assumption, 

F^ (P, «) = Fl{v) + /(«)[^°(P) - Fl{v)\ (3) 

where Fl{v) = Ff\'p, a = 1) = 1 + k - k/(1 + fi'^^p/n) 
andi^O(p) = F^"'(p,a = 0) = l-hK-Av/(l-F(At°'^p-hc)/K). 
F"'' {p, a) interpolates between F° (p) and F^ (p) through 
/(a) = (1 — a^)^/(l + 0"^ +ck^), which is chosen to guar- 
antee for the functional the second order gradient ex- 
pansion, good exchange jellium surface energies, and the 
Hartree-Fock exchange energy of the 12-electron hydro- 
genic density (ls^2s^2p^3s^) with the nuclear charge Z = 
1 [23] {Z — Im the following discussions if not mentioned 
otherwise) . See Ref. [2J for the derivatives of F™^ (p, a) 
with respect to p and a which are needed for the self- 
consistent implementation. For a slowly varying density 
where a « 1, the second term of the left hand side of Eq.[2] 
is negligible and of order V^ . F^"* {p, a) then reduces to 
F,}{p) and recovers the second order gradient expansion 
with the first principle coefficient /i*^^ = 10/81 f25'| as 
in PBEsol ,6j. c = 0.28771 and n = 0.29 are two pa- 
rameters fixed by the exchange energies of the hydrogen 
atom where a = and the 12-electron hydrogenic density 
which is used to guide the functional from a — Oto a ^ 1. 
They also deliver excellent exchange energies for other 
hydrogenic densities (see Ref. |24I ). No point in the (c, 
k) parameter space could be found without violating the 
loose Lieb-Oxford bound [H (ie., k sC k^^ = 0.804 0) 
if we use a spin-unpolarized hydrogenic density with elec- 
tron number less than 12. The obtained k — 0.29 sug- 
gests a tight Lieb-Oxford bound and leads to a very 
fiat F^{p) as shown in Fig. 1(b) The resulting small 

derivative at nuclei, — ^ d^'"~ ls=o.376 — 0.22, comes 
close to satisfying the finiteness constraint on the ex- 
change potential at nuclei. Compared to revTPSS, this 
much simpler form doesn't have the order-of-limits prob- 
lem while recovering regions of small s around a ~ I oi 
revTPSS as shown in Fig. 1(a) In the following discus- 



sion, we combine this exchange functional with the vari- 
ant of the PBE correlation (denoted as vPBEc), where 
/3(rs) = 0.066725(1 + 0.1r^)/(l 4- 0.1778rs) as used in 
revTPSS Q. Although the vPBEc is not one-electron 
self-correlation free, its error is small (about 0.006 Ha for 
the H atom), which is usually largely cancelled out of 
atomization energies involving H. It has an accurate cor- 
relation energy for a two-electron ion of nuclear charge 
Z ^ oo, which is -0.0479 Q, better than -0.0510 from 
TPSS and -0.0527 from revTPSS in comparison with the 
exact value -0.0467 [2J|. And it also has accurate corre- 
lation jellium surface energies Q. The resulting MGGA 
respects a tight Lieb-Oxford bound with Fxc ^ 2.137, 
while the loose Lieb-Oxford bound is 2.273 [il,[ii|. 

From the experience of GGAs [a,[a], we know the faster 
an enhancement factor grows with s, the more preference 





(b) 



FIG. 1: Exchange enhancement factors vs. a (a), and s (b). 



of the functional towards less compact systems. Since 
the enhancement factor of the present exchange func- 
tional as shown in Fig. 1(b) is largely depressed for a 



large range of s compared to that of revTPSS, it then 
seems to be a reasonable guess from the experience that 
the present exchange functional with the vPBEc would 
give too small lattice constants for solids and too high 
atomization energies for molecules. However, our results 
(given later) show a different scenario that our MGGA 
performs equally well for atoms, molecules, surfaces, and 
solids. The seeming contradiction between the prefor- 
mance and Fig. 1(b) is well resolved by resorting to the 
a-dependence shown in Fig. l(a)[ which is much less un- 
derstood. In previous constructions and analysis of MG- 
GAs [i, H, [Tj], it is emphasized that a enables MGGAs 
to distinguish the single-orbital regions from the orbital- 
overlap regions. Here, we further stress that the mono- 
tonically decreasing dependence of an enhancement fac- 
tor on a is qualitatively equivalent to the monotonically 
increasing s-dependence, which explains the seeming con- 
tradiction and will be rationalized by the following two 
observations. 

The first is on the changes in s and a distributions 
from the 10-electron hydrogenic density {ls^2s^2p^) to 
the 12-electron one {ls^2s^2p^3s^) ^, and on the cor- 
relation between the changes. Note these two densities 
are spherically symmetric. Fig. [5] shows the s and a 
distributions of these two densities. The shell structure 
can be easily recognized and roughly identified by a with 
a < 1 for shell regions and a > 1 for intershell regions. 
We can tolerate the confusion caused by this definition 
for the tail regions, where a could be or ^ 1 as shown 
in Fig. [21 since the tail regions are energetically less im- 
portant. When an electron is in the intershell region, and 




FIG. 2: The s and a distributions for the 10- and 12-electron 
hydrogenic densities. The line y=l is plotted to help roughly 
identify the shell (a < 1) and intershell (q > 1) regions. 



a is large even where s is sniall, the electron's exchange 
hole is probably not centered close to the electron, but 
is spread out over the smaller inner shell and the larger 
outer shell. This spreading of the hole will make the 
exchange energy density in the intershell region less neg- 
ative than it would be for a slowly-varying or uniform 
density. One can imagine i^^; < 1, as can happen in the 
present exchange functional for small s and large a that 



is shown in Fig. 1(a) 



When two 3s hydrogenic electrons are added to the 10- 
electron hydrogenic density, part of the outermost shell 
region (a < 1) changes into the intershell region (a > 1), 
which is associated with a decrease of s. Although it's 
not certain that an increase of a is always associated 
with a decrease of s during the formation of the inter- 
shell of an atom, it's very likely for the intershell region 
between the outermost core and the valence of an atom 
within a solid, an important region for determining the 
lattice constants of solids [2y, [27| (as will be discussed in 
the second observation). Then, the presumable correla- 
tion between s and a in the intershell regions of a solid 
suggests that monotonically decreasing a-dependence of 
an enhancement factor has the qualitatively same effect 
as monotonically increasing s-dependence does for these 
regions. 



TABLE L Error statistics of lattice constants (A) of the SL20 
solids J28| ] and atomization energies (kcal/mol) of the AE6 set 



See the text for the definitions of F^ 



F* 



and F" 





LSDA PBEsol 


^Omt 


^*mt 


F^' 


ME 
SL20 MAE 


-0.081 -0.012 
0.081 0.036 


-0.079 
0.079 


-0.126 
0.126 


0.016 
0.023 


ME 
AE6 MAE 


77.4 35.9 
77.4 35.9 


55.1 
55.1 


22.0 
22.6 


0.6 
5.5 



The second observation concerns the variations of the 
lattice constants of the set of 20 solids (SL20) [2^ and 
the atomization energies of the AE6 molecule set [2^ 
in response to changes of the a dependence of the en- 
hancement factor. In order to show the dramatic effect of 



the a dependence on the structural and thermochemichal 
properties and to deduce its origin, f{a) is compared 
to two variants. The first variant is /°(a) — 0, which 
is independent of a and actually results in the PBEsol 
GGA but with different k and (3{rs)- The second one is 
f*{a) = /(a) * hia) with h{a) = 2/[e("-i)/'' + 1] - 1. 
h(a) is equal to 1 for a < 1 and to -1 for a > 1, re- 
spectively, if 7 ^- 0. Here, we choose 7 = 0.1 for nu- 
merical reasons. Compared to /(a), /*(«) flips the a 
dependence for a > 1 from monotonically decreasing to 
monotonically increasing, and therefore favors more the 
regions with a > 1. The choice of the demarcation point 
at a = 1 is natural in view of the first observation and 
because it helps to satisfy the second order gradient ex- 
pansion, /(a) and the two variants — whose curves as 
functions of a are given in Ref. |2J — result in three dif- 
ferent enhancement factors, F™*, ^i?'"', and f^'"*. 



TABLE II: Error statistics for various density functionals. 
The LSDA and PBE enthalpies of formation are from Ref. |y. 
The lattice constants for M06L are from Ref. 30i, where 18 
solids, mostly overlapping with the SL20 set, were tested, a: 
Ref. [33 



LSDA PBE 



M06L 



revTPSS present 



Exchange energies (Ha) of rare gas atoms 
ME 2.274 0.219 0.204 0.291 0.068 

MAE 2.274 0.219 0.210 0.293 0.111 

Atomization energies (kcal/mol) of the AE6 molecules 
ME 77.4 12.4 3.2 3.3 0.6 

MAE 77.4 15.5 4.2 5.9 5.5 

Dissociation energies (kcal/mol) of the W6 water clusters 
ME 5.2 0.0 -0.2 -1.0 0.0 

MAE 5.2 0.3 0.4 1.0 0.1 

Enthalpies of formation (kcal/mol) of the G3 molecules 
ME -121.9 -21.7 -1.6 -3.6 -1.6 

MAE 121.9 22.2 5.2 4.8 8.3 

Jellium surface exchange energies (%) 
MRE 45.8 -20.9 -75.9 -1.0 -7.3 

MARE 45.8 20.9 75.9 2.2 8.0 

Jellium surface exchange-correlation energies (%) 
MRE -0.4 -3.1 24.5 2.6 -0.3 

MARE 0.4 3.1 24.5 2.6 1.6 

Lattice constants (A) of the SL20 solids 
ME -0.081 0.051 0.015" 0.015 0.016 

MAE 0.081 0.059 0.071" 0.033 0.023 



TableUshows the mean error (ME) and the mean abso- 
lute error (MAE) of the lattice constants of the SL20 set 
[28| . and the atomization energies of the AE6 set [20| 
from LSDA, PBEsol, and variants of F™'. The alle- 
viation, from LSDA to F°'"^^ and then to F^'^', of the 
overestimation in the atomization energies and of the 
underestimation in the lattice constants, suggests that 
the built-in monotonically increasing s-dependence and 
monotonically decreasing a-dependence in the enhance- 
ment factors reduce the preference of LSDA towards com- 
pact systems. However, from F"* to F^'"', where only 



the monotonicity of the a-dependence in the range of [1, 
oo] is flipped, the sohds in the SL20 set are drastically 
shrunk to such a surprising degree that the lattice con- 
stants of i^^'"' are significantly smaller than even those 
of LSDA. Since the important region in terms of deter- 
mining the lattice constant of a solid for a functional has 
been identified [2a, [23] to be the intershell region of the 
constituent atoms between the outermost core and the 
valence regions, the drastic shrinkage from F^"' to F*'"* 
is a strong indication that the sheh and intersheU re- 
gions are associated with q < 1 and a > 1, respectively. 
Compressing a sohd turns part of the outermost core 
and the valence regions (a < 1) into intershell regions 
(a > 1) between them, which F*'"' favors more than 
F^"' does. The absence of the monotonicahy decreasing 
a-dependence in F^?'"*, leading to the shrinkage of the 
solids, could be compensated by enhancing the mono- 
tonically increasing dependence on s, as PBEsol does by 
using K^'^ . This implies a decrease of s during the for- 
mation of the intershell regions and thus corroborates the 
correlation between s and a during the formation of the 
intershell region observed in Fig. [51 Therefore, in F"*, 
both the monotonically increasing s-dependence and the 
monotonically decreasing a-dependence have the effect 
of penalizing the formation of the intershell regions and 
enlarging the lattice constants. Similar deterioration is 
also found for the atomization energies of the AEG set for 
F*''^* compared to F^"*. Unlike for solids, F*'"' still sig- 
nificantly improves the atomization energies of the AE6 
set over F"'"', and thus LSDA, suggesting that the a- 
dependence in the range of [0, 1] has stronger influence 
in atoms and molecules than in solids. Remarkably, the 
monotonically decreasing a-dependence used in F"* sig- 
nificantly improves the overestimated atomization ener- 
gies of the AEG set of F°'"' to an excellent accuracy level 
with MAE of 5.4 kcal/mol. This implies that monotoni- 
cally decreasing a-dependence is in general able to make 
a functional favor less compact systems, as does mono- 
tonically increasing s-dependence. 

Now, let's turn to the results for atoms, molecules, sur- 
faces, and solids, which are summarized briefly in Table [U 
in terms of the ME and MAE, or their relative analogs 
MRE and MARE. See Ref. M for full details of Table [fil 
Table [ni shows that the use of the exact exchange energy 
of the 12-electron hydrogenic density guarantees excellent 
exchange energies for atoms, resulting in good atomiza- 
tion energies for this simple functional. In all categories 
shown in Table [Til the present functional outperforms 
the standard PBE GGA (PBE is shghtly better than 
the present functional for the cohesive energies of the 
20 SL20 sohds as shown in Ref.'2i). Within the MGGA 



level. The heavily parameterized MOGL [T3| predicts ex- 
cellent atomization energies and dissociation energies for 
the W6 set, at which it aims during the construction. 
However, it is significantly wrong for the jellium sur- 
face energies contributed from the exchange and correla- 
tion terms, separately or together. The too-large MOGL 
jellium surface exchange-correlation energies imply that 
metal bulks are overstabilized, consistent with the too- 
small lattice constants of main group simple metals, e.g., 
Na and Al \30\. The present functional and revTPSS are 
more balanced for different categories and therefore more 
robust. Compared to revTPSS, the present functional is 
better for the SL20 solids, comparable for the jellium 
surface exchange-correlation energies but worse for the 
exchange part alone, and worse for the G3 molecules. 
However, the present functional predicts the most accu- 
rate dissociation energies of the W6 water clusters among 
the functionals, implying a good description for the hy- 
drogen bond. 

In summary, we have for the first time stud- 
ied the effect and importance of the dependence of 
computationally-efficient semilocal MGGAs on the ki- 
netic energy density through the dimensionless inhomo- 
geneity parameter a, and presented a new MGGA ex- 
change functional that disentangles a from the reduced 
density gradient s by the means of separability assump- 
tion. By varying the a— dependence in the exchange 
functional, we showed that the formation of the inter- 
shell region between the outermost core and the valence 
of an atom within a solid is associated with an increase 
of a and a decrease of s, suggesting that monotoni- 
cally decreasing a-dependence of an enhancement fac- 
tor is qualitatively equivalent to monotonically increas- 
ing s-dependence for these intershell regions. This has 
a significant impact on the construction of MGGAs and 
the MGGA-based nonlocal approximations, as exempli- 
fied by the present MGGA — which is overall comparable 
in performance, but quite different and much simpler in 
form, compared to the sophisticated revTPSS MGGA, 
and thus demonstrates the fiexibility and the rich struc- 
ture of MGGA brought by the extra ingredient of the 
kinetic energy density. 
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